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Abstract 

In this report, the implementation of a program for the simulation of the sky 
map, in the gamma rays energy range, observed by the GLAST experiment 
will be described. The program generates a list of photons and images of 
the galactic and extragalactic background and of the sources, in a selected 
energy range and in a given region of the sky. 

1 Introduction 

The simulation program is organized as follows: the galactic background map can 
be generated using the GALPROP program P or the model of the diffuse galac¬ 
tic background [21 obtained using observations from EGRET. The extragalactic 
contribution is given by a constant value in a fixed energy range. The gamma 
emission of the sources is parameterised using a standard power law and sources 
from the Third Egret Catalog jSj and faint sources generated following the Stecker 
and Salamon model @1 are considered. All the contributions (background and 
sources) are integrated in the given energy range {[Emin, Emax]), in a fixed re¬ 
gion {Ab X A/) of latitude and longitude and convoluted with the instrumental 
point spread function (PSE), effective area (SA) and energy resolution (ED) of 
the LAT instrument [Hj. Einally the total intensity is multiplied by the exposure 
time. 

2 Background Description 

Eor the Gamma Ray Galactic Background either GALPROP program or the 
EGRET model can be used. GALPROP is a simulation program based on a 
model which reproduces many kinds of observational data related to cosmic-ray 
origin and propagation. This model provides a good basis for studies of galactic 
gamma-ray background taking into account the contributions from the most im¬ 
portant emission mechanisms. The other possibility is to use the EGRET map 
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Figure 1 : Contribution to the photon flux in the region centred on 3C279 after 
one precession period in the energy range 0.1 ^30 GeV. Left: sources only. Right: 
background plus sources. The colour bar indicates the number of counts. 


measured in the energy range between 0.1 and 30 GeV. The diffuse extra galactic 
background can be described using the differential flux of photons ^ : 

^ = 7.3 • 10“® • 0.4512-1 • (photons ■ cm-‘^ ■ s"! • GeC"! • sfi) (1) 
dE 

The total contribution of the extra galactic background is obtained by integrating 
Equation Q between and Emax- 

3 Detector Simulation 

The Effective Area (SA), Point Spread Function (PSF) and Energy Resolution 
(ED) of the detector are energy dependent, the dependence of the effective area 
on the energy E and on the photon incidence angle 0 can be parameterised as: 
SAiO^E) = SA1{E) • SA2{0). For a position at 7 ^ 0 there is, as a first approx¬ 
imation, also a dependence of the effective area on the angle. The PSF can be 
described using two different functions: in the simplest case the PSF (for normal 
incident gammas) is assumed to be Gaussian. The radius containing the 68 % of 
the photons is taken as the RMS of the PSF distribution. As a first approximation 
we assume that cFpsF is independent on the incidence angle of the photons, but 
this hypothesis does not take into account that most of the exposure comes from 
fairly far off axis, where the PSF is broader. A better approximation is to consider 
a superposition of two Gaussians. Even if this choice underestimates the broad 
tails, it takes into account both the narrow core and the broad component of the 
PSF. The energy spread function is also assumed to be Gaussian and, as a first 
approximation, independent on the inclination. 
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4 Orbit and Exposure Time Simulation 

The duration of the exposure can be either fixed, or calculated using an orbit 
simulator. The orbit is simulated with a geometrical approach. Orbital parame¬ 
ters, satellite positions and zenith pointing directions are evaluated in steps of 30 
s 2°) and assuming that within each step the orbit parameters are constant. 
The total exposure is calculated adding the contributions over all the steps. 

If the total observation time t is not fixed, the exposure £ is computed for each 
point of the sky in a given region A{Ab x A/) simulating the orbit of the space¬ 
craft (SC). The exposure will be zero if the SC is in the South Atlantic Anomaly 
(SAA) or if the angle between the zenith direction and the source is greater than 
105°. Outside these regions it will be the product of the exposure time by the 
effective area SA at energy E. If the time interval t is fixed, the exposure is 
£{E,0) =SA{E,0)^t. 

5 Convolution with PSF, SA and ED 

The differential flux C(Po) (photons per unit of area, time and solid angle), is given 
by the sum of all the contributions: the constant extragalactic background, the 
galactic one and the sources. The output is a map of the sky that at this stage 
does not yet include detector resolution effects, in order to have the real image 
observed by the LAT experiment, we have to convolute our result with the SA, 
PSF and ED. We define a distribution of the photons as function of the energy 
given by D{E) = £{E,0) • I{E)dO. The flux I{E) is calculated from: 

I{E) = /o • E-^ {photons • cm-2 • 5“^ • GeE"^) (2) 

where / is the flux intensity , E is the photon energy, Iq is a constant taken from 
the Egret Catalog and a is the spectral index, assumed to be equal to 2.1 for 
the background. The total number N' of photons is obtained by integrating the 
function D{E) between Emin and Emax'- the generated number of photons is then 
selected assuming a Poissonian statistics. The N photons detected are distributed 
inside the sky region assigning to them a random energy Etme according to the 
D{E) distribution and an inclination angle distributed according to the function 
£{Etrue^ ^)- Then the energy measured by the detector, {Emeas)i is obtained using 
the function ED {Etme: Emeas) and the angular distance, p, from the origin point 
Po = {h: bo) consistent with the PSE{p, Etme)- The arrival time of the photons is 
generated taking into account the Poissonian nature of the process and the LAT 
dead time (assumed to be 100/is). 

As a first approximation the sources are assumed to be point-like, located at an 
infinite distance and emitting an intensity of photons described by Equation [2l 
Faint sources can be considered generating their flux P (for E > 100 MeV) in the 
range between photons • cm~‘^ • s~^ (GLAST lower limit) and photons • 
cm~‘^ • s~^ (EGRET lower limit) according to: 


N = No’E^ 


( 3 ) 
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Figure 2: Total number of photons, in a logarithmic scale, from sources (right) 
and background plus sources (left) after one precession period in the energy range 


0.1^30 GeV. 


where N is the number of sources with flux F, and A^o and a are extrapolated 
using the predicted distribution from Stecker and Salamon. In the case of faint 
sources the spectral index in Equation (21 is generated according to a Gaussian 
distribution centred in 2.1 with an RMS equal to 0.2 (classic range for Blazar 
spectral index). The total number of sources is obtained integrating Equation IHI 
over the flux range and weighting for the observed fraction of solid angle. The 
total number of photons from a source is calculated from Equation [2l and they 
are distributed in the map following the same procedure as described above for 
the background. In Figure Q the map with the contribution from the sources and 
background, is shown, while Figure O shows a simulation of the whole sky. 


6 Results and Conclusions 


The result of the simulation yields images and two tables containing the informa¬ 
tion required for the analysis. Parameters like the position of the sun and moon 
have been derived using astrophysical routines which calculates the position of 
the planets in the solar system knowing the Julian date; while the geomagnetic 
coordinates have been computed using the GEOPACK [Jj code. The output is 
produced at different stages of the simulation, where the time step is given as 
input parameter. 

A simulation program to calculate the flux of photons observed by the GLAST 
experiment has been implemeted. Contribution from the galactic, extra-galactic 
background and sources has been included. The effect of the detector is included 
by convoluting the obtained fluxes with the Point Spread Function (PSF), Effec¬ 
tive Area (SA) and Energy Resolution (ED). The time is given either as parameter 
or it is calculated from the simulation of the orbit. 








C. Cecchi A Fast Simulator for the GLAST experiment 


5 


References 

[1] A.W. Strong, LV. Moskalenko, New Developments in the GALPROP CR 
Propagation Model astro-ph/0106504, (2001) 

[2] S.D. Hunter et dl.^EGRET Observations of the Diffuse Gamma Ray Emission 
from the Galaetie Plane Astrophys. J., (1997) 

[3] R.C. Hartman et ak. The Third Egret Gatalog of High Energy Gamma Ray 
Sourees Astrophys.J.SuppL, (1999) 

[4] F.W. Stecker, M.H. Salamon, The Gamma Ray Baekground from Blazars: a 
New Look Phys. Rew. Lett, (1996) 

[5] http://glast.stanford.edu 

[6] P. Sreekumar, EGRET Observations of the Extragalaetie Gamma Ray Emis¬ 
sion Astrophys. J., (1998) 

[7] N.A. Tsyganenko, Tsyganenko Magnetie Eield Model and GEOPAGK s/w, 
http: / / nssdc.gsfc.nasa.gov/space/model/magnetos/tsygan.html 



